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Abstract 

We study the minimum volume ellipsoid estimator associates to a 
cloud of points in phase space. Using as a natural measure of uncer- 
tainty the symplectic capacity of the covariance ellipsoid we find that 
classical uncertainties obey relations similar to those found in non- 
standard quantum mechanics. 
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1 Introduction 

Contrarily to what is often believed the Heisenberg uncertainty principle 

ApjAxj > \h (1) 

is not a statement about the accuracy of our measurement instruments; 
its derivation assumes on the contrary perfect instruments. The correct 
interpretation of Heisenberg's inequalities is the following (see e.g. Peres 
[15] . p. 93): if the same preparation procedure is repeated a large number of 



* E-mail: maurice.de. gosson@univie.ac. at 



1 



times, and is followed by either by a measurement of Xj, or by a measurement 
of pj, the results obtained have standard deviations Axj and Apj satisfying 
(fT|). The same interpretation is of course true for the stronger Robertson- 
Schrodinger [T71 [21] inequalities 

{A Pj ) 2 {A Xj f > A{x vPj ) 2 + \h 2 (2) 

to which ([1]) reduce if one neglects the covariances A(xj,pj) 2 ; they are com- 
plemented by the trivial inequalities 

ApjAx k > if j ^ k, ApjApk > , AxjAxk > (3) 

which might be violated in a nonstandard form of quantum mechanics, 
noncommutative mechanics (NCQM, see Dias et al. [3] and the references 
therein), where the second and third inequalities ([3]) are be replaced with 

Ap 2 Ap 2 >A( Pj , Pk ) 2 + \9% , Ax 2 Ax 2 k >A( Xj ,x k ) 2 + \r, 2 k (4) 

where 6 2 k = 9 2 kj and rfj k = r] 2 kj - 

In classical statistical mechanics the situation is somewhat different: due 
to the inherent inaccuracy of the measurement apparatus there are uncer- 
tainties for all pairs of variables, conjugate or not. We are going to show in 
this article that despite the different nature of quantum and classical uncer- 
tainties, the latter can be described in a similar way, leading to inequalities 
which are formally identical to those of NCQM, namely 

Ax 2 Ax 2 k > A( Xj , x k ) 2 + a 2 k (5) 
Ap 2 Apl>A( Pj ,p k ) 2 + c 2 k (6) 
Ax 2 Ap 2 k > A{ Xj ,p k ) 2 + b 2 k (7) 

with a 2 k = a k j, b 2 k = 6|-, c 2 k = c\-. We have actually already shown in a 
recent paper [S] (also see de Gosson and Luef [7]) that inequalities of the type 
(PQ)-([2]) are by no means characteristic of quantum mechanics, and that there 
are formally similar statements in classical statistical mechanics; to sustain 
our claim we used tools from robust multivariate statistics (the Minimum 
Volume Ellipsoid method, reviewed in Section [2]) together with a topological 
device, the notion of symplectic capacity, which we use as a natural device 
for measuring uncertainty. The symplectic capacity of a closed phase space 
set is the equatorial area of the largest ball that can be sent in this set using 
only symplectomorphisms (=canonical transformations). The symplectic 
capacity of the covariance ellipsoid is not related to its volume (it always has 
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the dimension of an area), but rather to that of a classical action. We note 
that symplectic capacities have been used by Scheeres and his collaborators 
[TTJ[2U] to study satellite guidance and constraints on spacecraft trajectories. 

The aim of this paper is to extend the results in [5] and to show that 
these NCQM uncertainties also appear quite naturally in classical statistical 
mechanics. 

Notation 1 We identify R n x R n with R 2n and write (x,p) = z if x <E W 1 
and p € M. n . We will view x,p and z as column vectors in all matrix 
computations. A symplectic form on IR 2 ™ is a non- degenerate bilinear an- 
tisymmetric form on that space. Let a be a symplectic form on M. 2n ; the 
corresponding symplectic group is denoted Sp(2n, a) [it is the group of all 
linear automorphisms S ofM? n such that a(Sz,Sz') = a(z,z') for all vec- 
tors z,z' in IR 2 ™/. We denote by a the standard symplectic form on M? n : 
a(z,z') = (z') T Jz if z = (x,p), z' = (x',p'); J is the standard symplectic 



2 The Minimum Volume Ellipsoid Method 

Let us consider a system of point-like particles; we assume that there are 
n degrees of freedom and label the position coordinates and momenta x = 
(xi, x n ), p = (p\, ...,p n ); we will also use the collective notation z = (x,p). 
Assume now that we perform position and momentum measurements on a 
large number K of identical copies of that system; we get a cloud S = 
{zi, Z2, zn}, N = nK, of points in M 2n . An efficient method for studying 
that cloud consists in using the minimum volume ellipsoid (MVE) method 
for the multivariate location and scatter (Rousseeuw [IB])- Geometrically 
speaking this method is an application of the John-Lowner theorem (see for 
instance Ball [2]): consider a subset {z^ , Zi 2 , Zi k } of >S; we will assume that 
the points z%. are in general position, i.e. that they do not remain in some 
hyperplane of M 2n . This condition is sufficient and necessary for any ellipsoid 
containing these points to have positive volume. The points %i x , Z{ 2 , Z{ h 
determine a polyhedron in R 2n ; we denote by fCk the convex hull of that 
polyhedron (it is the smallest subset of M 2n containing {2^ , Zi 2 , Zi k }). The 
John-Lowner theorem ensures us that there exists a unique ellipsoid in 
M 2n containing /C^ and having minimum volume among all the ellipsoids 
having this property. Repeating this process for all subsets of the cloud 
«S having k elements in general position we get a family of ellipsoids; by 
definition the MVE is the one with the smallest volume. 
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This method is practically implemented as follows (see e.g. Van Aelst 
and Rousseeuw [22], Hubert et al. |12|): choose an integer k between \N/2] + 
1 and ./V ([N/2] the integer part of N/2); this constant k determines the 
robustness of the resulting estimator; a common choice is 



k 



N + 2n + 1 



(8) 



By definition, the location and scatter estimators minimize the determinant 
of the matrices M subject to the condition 

# {j : (Zj - z) T M-\ Zj -z)<m 2 }>k (9) 

where the minimization is over all z € M? n and all positive definite symmetric 
matrices M of size 2n. Here to is a fixed constant, chosen so that the MVE 
estimator is a consistent estimator for of the covariance matrix for data 
coming from a multivariate normal distribution, that is 



m = yxL, a , a = k / N 

where xin Q is a chi-square distribution with 2n degrees of freedom (see 
Lopuhaa and Rousseeuw |14|). One the pair (M,z) is determined, the min- 
imum volume ellipsoid (MVE) is the set of all z in M? n such that 

[z - z) T M~ l (z - z) < m 2 . (10) 

The next step consists in associating to the MVE J a covariance ma- 
trix. For this one has to choose an adequate value mo for to; denoting the 
corresponding matrix M by S the MVE is the ellipsoid 

C : {z-zf^iz-z) < ml (11) 

and S is then precisely the covariance matrix. 

We will write the covariance matrix in the form 



A(x,x) A(x,p) 
A(p,x) A(p,p) 



(12) 



where A(x,x) = ((A(xj, Xj)) lKi - <n and so on. We will write, as is custom- 
ary, A(xj,Xj) = Axf and A(pi,pi) = Apf; we have A(xi,xj) = A(xj,Xi), 
A(pi,pj) = A(j)j,pi), A(xi,pj) = A(pj,Xi) hence the covariance matrix is 
symmetric. 
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3 A Condition on E 



Let A = (a,jk)i<j : k<n and C = {cjk)i<j,k<n be two real antisymmetric matri- 
ces, and B = {bjk)i<j,k<n a rea l symmetric matrix. To A, B, C we associate 
the In x 2n antisymmetric matrix 







A B 
-B C 



which is the most general form an antisymmetric of size 2n can have. We will 
always assume that the matrix Q is in addition invertible. This condition 
implies that the bilinear form u on M. 2n defined by uj(z, z') = —(z') T Q,~ 1 z 
is a symplectic form on IR 2n , and we denote by (M. 2n ,uj) the corresponding 
symplectic phase space. Notice that when A = C = and B = I we have 
I> 



Q = J = I j. q I , the standard symplectic matrix. The general case can 
actually be reduced to the standard case: 

Proposition 2 There exist linear automorphisms F of M? n such that fl = 
F T JF; equivalently F is a linear symplectomorphism (IR 2n , a) — > (R. 2n ,uj), 
that is we have uj(Fz, Fz') = a(z,z') for all z and z' . 

Proof. Let £ CT = {ef , <}U{/f , /£} and B» = {ef, <}U{/f , /«} 
be symplectic bases of (M 2n ,cr) and (M. 2n ,uj) respectively (i.e. a;(e^,e^) = 
= Oi <jJ (fj > ^ e k) = ^j'fc an d similar relations for the <r(ej,e£), etc.). 
The automorphism F defined by F{e") = e^ , F{fJ) = fj 3 is a symplecto- 
morphism (M 2n , cr) — > (R 2n ,u). m 

We remark that this result (which can also be proved using the properties 
of the Pfaffian Pf(Q)) is just a restatement of the linear (and global) version 
of Darboux's theorem [9] on the local equivalence of all symplectic manifolds 
with same dimension. 

Let now E be the covariance matrix as defined above, and consider the 
matrix 

V4-;o- { A (x,x) + iA A(x,p) + iB\ 
+ \A(p,x)-iB A(p,p)+iCj- l j 

We observe that E + i£l is Hermitian since E is symmetric and (iCl)* = iVt 
since is real antisymmetric. The eigenvalues of E + iQ, are thus real. From 
now on we will assume that these eigenvalue are nonnegative, that is E + iQ, 
is semi-definite positive, which we write 

£ + *fi>0. (14) 
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(One can actually show that this condition automatically implies that S > 
0). Notice that we have in particular 

A(x,x) + iA>0 , A(p,p) + iC > 0. (15) 

When n = 1 the covariance matrix is just 

Ax 2 A(x,p) N 
A(p,x) Ap 2 

and the antisymmetric matrices and ./V are zero so that H = aJ 

1 The condition £ + iQ, = £ + iaJ > is in this case equiva- 



a 
lent to 

Ax 2 + za A(x,p) \ > p 

A(p,x) Ap 2 — ia J ~ 

which is in turn equivalent to the single inequality 

Ax 2 Ap 2 > A(x,p) 2 + a 2 . 

This is of course formally the Robertson-Schrodinger inequality ([2]); in par- 
ticular we have the Heisenberg-type inequality AxAp > a. 

Let us extend the study to higher dimensions. When n = 2 the matrices 
A, B, and C are of the type 

a ( a \ n fb d\ „ ( c 
A= U 0j^=U eJ' C= U 

so that £ + zfi is the 4x4 matrix 

^ Axf A(xi,x 2 ) + ia A(xi,pi) + ib A(xi,p 2 ) + ict\ 

A(x2,xi) — ia Ax 2 . A(x 2 ,pi)+id A(x 2 ,p 2 ) + ie 

A(pi,a?i) — ib A(pi,x 2 )-id Ap\ A(pi,p 2 ) + ic 

\A(p 2 ,xi)-id A(p 2 ,x 2 )-ie A(p 2 ,p 1 )-ic Ap\ ) 

Recalling Sylvester's criterion [1CH 03] which says that a Hermitian matrix is 
positive semidefinite if an only if all of its principal minors are nonnegative, 
the condition £ + i$7 > implies that the principal minors of order two of 
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S + iCl must be > 0, we immediately get the inequalities 

AxjAx 2 . > A(xi,x 2 ) 2 + a 2 
AxjAp 2 > A(xi,pi) 2 + b 2 
ApfAp 2 2 > A( Pl , P2 ) 2 + c 2 
AxjAp 2 > A( Xl , P2 ) 2 + d 2 
AxjApl > A(x 2 ,pi) 2 + d 2 
Ax\Ap\ > A(x 2 ,p 2 ) 2 + e 2 . 

The same argument shows that more generally: 

Proposition 3 Let n > 2. If the covariance matrix £ satisfies the condition 
£ + iO, > then the following uncertainty relations hold: 

Ax 2 Ax 2 k > A( Xj ,x k ) 2 + a 2 k (17) 
Ap 2 Ap 2 >A{p 3 ,p k ) 2 + c 2 k (18) 
Ax 2 Ap 2 k >A{x v p k ) 2 + b 2 k . (19) 

In particular, if £1 = e J, e > 0, these conditions reduce to the Robertson- 
Schrodinger inequalities 

AxjAp 2 k > Aix^pkf + e 2 . 

A warning: the group of inequalities f)lT|) — (|19|) is not equivalent to the 
condition S + iVL > as soon as n > 1 . Here is a simple counterexample in 
the case n = 2 and Q = J: choose n = 1 and 



E = 



/ 1 


-1 







-1 


1 














1 













V 



This matrix is positive definite, and the inequalities above hold trivially 
(they reduce to equalities); the matrix S + iJ is however indefinite (its 
determinant is —1). 

4 The Symplectic Capacity of an Ellipsoid 

For the basics of symplectic geometry we are going to use we refer to Arnol'd 
PQ, de Gosson [1] (Arnol'd uses the term "canonical transformation" for 
symplectomorphism; this terminology is usual in Physics). 
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The condition £ + id > can be restated in terms of the symplectic 
capacity of the minimum volume ellipsoid Co given by (jlip . Let us recall 
the definition of the notion of symplectic capacity (of which we have given 
a detailed discussion in de Gosson and Luef [7j). Let us call symplectic 
manifold a submanifold U of M? n (possibly with boundary) equipped with 
a symplectic form a. A symplectic capacity associates to every symplectic 
manifold (U, a) a number c(U, a) > or +00; this correspondence must 
satisfy the following axioms [HI [16] : 

1. Monotonicity: If <3? : (U,a) — > (Y,/3) is a symplectic embedding, 
i.e. a diffeomorphism satisfying /3(<3?(z), $(z') = a(z, z') we must have 
c(U,a) <c(V,(3); 

2. Conformality: For every real A/Owe have c(U, \a) = X 2 c(U, a); 

3. Normalization: c(B(R),a) = c(Zj(R),a) = ttR 2 ; here a is the stan- 
dard symplectic form, B(R) the ball \x\ 2 + \p\ 2 < R 2 and Zj(R) the 
cylinder x 2 + p 2 < R 2 . 

In general symplectic capacities are not related to the notion of volume; 
this is already clear from the normalization condition c(Zj(R),a) = ttR 2 
which shows that as soon as n > 1 the symplectic capacity of a region with 
infinite volume can be finite; in fact the conformality axiom shows that 
symplectic capacities behave under dilations. Also note that sym- 

plectic capacities are extensive quantities, i.e. they do not depend directly 
on dimension, as volume does. 

The monotonicity axiom implies that if there exists a symplectomor- 
phism (=symplectic diffeomorphism) <1> : (U, a) — > (V, /3) such that <$(£7) = 
V then 

c(U,a)=c($(y),P)=c(V,p) (20) 

hence symplectic capacities are symplectic invariants. The basic example 
of a symplectic capacity is Gromov's widtlJ3 cqr- It is defined as follows: 
let i?GR be the "symplectic radius" of U: it is the supremum of all radii 
of balls B(R) that can be embedded in U using symplectomorphisms of 
(M? n ,a). By definition the Gromov width of U is then cqr(^, cr) = vri?Q R 
(with cgr(U, a) = 00 if Rqb. = 00). The fact that the normalization axiom is 
satisfied follows from Gromov's non-squeezing theorem [8j, a deep property 
of symplectic topology which says that the ball B(R) cannot be embedded 
inside a cylinder Zj(r) with radius r < R. If U is a compact, connected 

1 Sometimes also called symplectic area. 
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and simply connected domain in the symplectic plane (R 2 , a) = (M 2 , — det) 
then cqr(U, a) is the area of U. There exist infinitely many symplectic 
capacities (and cgr is the smallest of all) but they all agree on ellipsoids 
(see e.g. [Sl[7])- We are going to give an explicit formula below, but let us 
first introduce the following notation and terminology. Let M be a positive 
definite 2n x In real matrix and consider the product JM. Its eigenvalues 
are those of the antisymmetric matrix M 1 ' 2 JM 1 / 2 and are thus of the type 
±i\ a j(M), j = 1, ...,n, with X a j(M) > 0. The numbers A CT j(M) are called 
the a- eigenvalues of M; up to a simultaneous reordering of the variables Xj 
and pj one can always assume that A CTi i(M) > • • • > A<j jTl (M). The ordered 
set 

Spec ff (M) = (A (T)1 (M),...,A (Jin (M)) 

is then called the a- spectrum of M (when it is understood that it is the 
standard symplectic structure which is used one speaks about symplectic 
eigenvalues and symplectic spectrum). One proves the following properties: 

M<M'=^ \a,j(M) < Xa,j(M') ,j = l, ...,n (21) 

where M < M' means that M' — M is semi-definite positive, and 

(A (T ,i(M- 1 ),...,A CT , n (M- 1 )) = (A ffi „(M)- 1 ,...,A CT ,i(M)- 1 ) (22) 

(see de Gosson [1], §8.3). Moreover there exists a symplectic matrix S € 
Sp(2n, a) such that 

5 T M5=( A M A ° M ) (23) 

where A M is the diagonal matrix diag(A CT) i(M) > ••• > A t7 , n (M)) (Williamson's 
diagonal form, see [H HI [23] ) . Formula (f23|) implies the following: let Ai 
be the phase space ellipsoid by inequality z T Mz < 1. The the inverse image 
S~ 1 (A4) has the normal form 

£>,,-(jio(*? +*?)<!• ( 24 ) 

One can of course replace the matrix J above by the antisymmetric non- 
degenerate matrix Q = F T JF (F is defined as in Proposition [2]) . Consider- 
ing as above the symplectic form oj(z, z') = — {z') T £l~ l z the w-spectrum of 
M is the decreasing sequence 

SpecjM) = (A Wil (M),...,A w , n (M)) 
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of positive numbers such that the ±iX U} j(M) (X UJ j(M) > 0) are the eigenval- 
ues of of fiM; the properties (|2ip and ()22[) hold mutatis mutandis, replacing 
the subscript a with u. We have: 

SpecjM) = Spec a (FMF T ). (25) 

To prove this it is sufficient to show that QM and J(FMF T ) have the same 
eigenvalues. Suppose OM = F T JFMz = Xz for some complex number 
A and z 7^ 0. This is equivalent to JFMz = X(F T )~ 1 z and hence to 
JFMF T [(F T ) _1 z] = A [{F T y l z\ which proves our claim. 

Proposition 4 Xei M be definite positive and consider the ellipsoids A4 + = 
{z : z T Mz < 1} and M~ = {z : z T M^ 1 z < 1} in R 2n . For every symplectic 
capacity c we have 

c(M + ,a)=Tr/X a>1 (M) , c(M-,o-)=AirX CT)n (M). (26) 

Proof. See for instance [9] for a proof of the first formula (|26p . The 
second formula follows from the first in view of the equality (I22p . ■ 

The first formula (|26p shows that the symplectic capacity of A4 is the 
area of the intersection of that ellipsoid with the x\,p% plane once it has been 
put in normal form (|24p . This again shows that the symplectic capacity is 
related to an area, and not to volume. In fact, using the invariance of the 
action form one can restate the first formula (I26p in the following way: 

c{M,a) = j> pdx = | j> pdx — xdp (27) 

•'7min •> Train 

where 7 is the shortest Hamiltonian orbit carried by the surface of the ellip- 
soid M (see de Gosson and Luef [7]); the symplectic capacity of an ellipsoid 
is thus explicitly expressed in terms of a dynamical action. 

5 Application to the Uncertainty Principle 

Let us now prove the main result of this article. Recall (formula (jlip ) that 
the covariance matrix is the set 

C = {z : {z-z) T Y,- l (z- z) <m 2 Q } (28) 

for some suitable choice of the real number ttiq. 
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Theorem 5 The condition Y, + iVt > for the covariance matrix £ is equiv- 
alent (for every symplectic capacity c) to the condition 

c(C, u) = c((F T )~ 1 C, a) > 7rm^A„(S) (29) 

where A n (E) is the largest a- eigenvalue of X. // \29\) is satisfied then we 
have 

Ax 2 Ax 2 k > A( Xj ,x k ) 2 + a 2 jk (30) 
Ap 2 Ap\ > A(pj,p k ) 2 + c 2 k (31) 
Ax 2 Ap 2 >A{x 3 ,p k ) 2 + h 2 k . (32) 

Proof. Since translations are symplectomorphisms in any symplectic struc- 
ture, the ellipsoid C has the same symplectic capacity as the centered ellip- 
soid z T Ti~ 1 z < tuq. The result now follows from Proposition [5] above with 
M = m 1 £- 1 . ■ 

6 Discussion 

The notion of symplectic capacity appears as a device allowing to measure 
in a geometrical way the size of the MVE in a new way. The inequalities 

Ax 2 Ax 2 k > A( Xj ,x k ) 2 + a 2 k 
Ap)Ap\ > A{ P] ,p k ) 2 + c 2 k 
Ax 2 Ap 2 k > A{ Xj , Pk ) 2 + b%. 

are not equivalent to the condition c(C, uj) > 7t?tt.qA^ but are implied by it. 
Therefore, c(C,lo) > 7rm,QA n (X) can be viewed as a stronger -but natural- 
version of the uncertainty principle. Its usefulness might very well come 
from the fact that the condition c(C,u) > 7rmQA ri (E) is invariant under 
arbitrary symplectic transformations (linear or not). It is in particular pre- 
served under Hamiltonian time evolution, since Hamiltonian flows consist of 
symplectomorphisms . 

Acknowledgement 6 This work has been financed by the Austrian Re- 
search Agency FWF (Projektnummer P20^2-N13). 



11 



References 

[1] Arnol'd V I, 1989 Mathematical Methods of Classical Mechanics, Grad- 
uate Texts in Mathematics, 2nd edition, Springer- Verlag 

[2] Ball K M, 1992 Ellipsoids of maximal volume in convex bodies, Geom. 
Dedicata 41(2), 241-250 

[3] Dias, N C, de Gosson, M, Luef, F, and Prata, J N, 2010, A deformation 
quantization theory for noncommutative quantum mechanics, J. Math. 
Phys. 51 072101 

[4] de Gosson M, 2006 Symplectic Geometry and Quantum Mechanics, 
Birkhauser, Basel, series "Operator Theory: Advances and Applica- 
tions" (subseries: "Advances in Partial Differential Equations"), Vol. 
166 

[5] de Gosson M, 2009 The Symplectic Camel and the Uncertainty Princi- 
ple: The Tip of an Iceberg?, Found. Phys. 99 194-214 

[6] de Gosson M and Luef F, 2007 Quantum States and Hardy's Formula- 
tion of the Uncertainty Principle: a Symplectic Approach, Lett. Math. 
Phys. 80 69-82 

[7] de Gosson M and Luef F, 2009 Symplectic Capacities and the Geometry 
of Uncertainty: the Irruption of Symplectic Topology in Classical and 
Quantum Mechanics, Physics Reports 484 131-179 

[8] Gromov M, 1985 Pseudoholomorphic curves in symplectic manifolds, 
Inventiones Mathematica 82 307-347 

[9] Hofer H and Zehnder E, 1994 Symplectic Invariants and Hamiltonian 
Dynamics, Birkhauser Advanced texts (Basler Lehrbiicher), Birkhauser 
Verlag 

[10] Horn R A and Johnson C R, 1985 Matrix analysis, Cambridge Univer- 
sity Press 

[11] Hsiao F.-Y and Scheeres D J, 2007 Fundamental Constraints on Uncer- 
tainty Relations in Hamiltonian Systems, IEEE Transactions In Auto- 
matic Control 52(4) 686-691 

[12] Hubert M, Rousseeuw, P J, and Van Aelst S, 2008 High-Breakdown 
Robust Multivariate Methods, Statistical Science 23(1) 92-119 



12 



[13] Johnson R C (ed.), 1990 Matrix Theory and Applications, Amer. Math. 
Soc, R I 

[14] Lopuhaa H P and Rousseeuw P J, 1991 Breakdown points of afhne 
equivariant estimators of multivariate location and covariance matrices, 
The Annals of Statistics 19(1) 229-248 

[15] Peres A, 1993 Quantum Theory: Concepts and Methods, Kluwer Aca- 
demic Publishers 

[16] Polterovich L, 2001 The Geometry of the Group of Symplectic Diffeo- 
morphisms, Lectures in Mathematics, Birkhauser 

[17] Robertson H P, 1929 The uncertainty principle, Phys. Rev. 34 163-164 

[18] Rousseeuw P J, 1985 Multivariate estimation with high breakdown point. 
In Grossmann, W, Pflug, G, Vicenze, I & Wertz, (Eds.), Mathematical 
Statistics and Applications, Vol. B., Riedel Publishing, Dordrecht the 
Netherlands 283-297 

[19] Rousseeuw P J and Leroy A M, 1987 Robust regression and outlier 
detection, John Wiley & Sons, New York 

[20] Scheeres D J, Hsiao F.-Y, Park R S, Villac B F, and Maruskin J M, 2006 
Fundamental Limits on Spacecraft Orbit Uncertainty and Distribution 
Propagation, Journal of the Astronautical Sciences 54 505-523 

[21] Schrodinger E, 1930 Zum Heisenbergschen Unscharfeprinzip, Berliner 
Berichte 296-303 [English translation: Angelow, A., Batoni, M.C.: 
About Heisenberg Uncertainty Relation. Bulg. Journal of Physics, 26, 
nos.5/6, 193-203 (1999), and http://arxiv.org/abs/quant-ph/9903100l 

[22] Van Aelst S and Rousseeuw P J, 2009 Minimum volume ellipsoid, Wiley 
Interdisciplinary Reviews: Computational Statistics 1(1) 71-82 

[23] J. Williamson, 1936 On the algebraic problem concerning the normal 
forms of linear dynamical systems, Amer. J. of Math. 58 141-163. 



13 



